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The Landau gauge lattice QCD simulation and the gluon propagator 
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The gluon propagator in the Landau gauge lattice QCD simulation is measured. The data suggests the con- 
finement mechanism of the Gribov-Zwanziger theory. 



1. Introduction 

In 1978 Gribov showed that the fixing of the di- 
vergence of the gauge field in non-Abelian gauge 
theory does not fix its gauge M. The restriction of 
the gauge field such that its Faddeev-Popov de- 
terminant is positive (Gribov region) is necessary 
yet insufficient for excluding the ambiguity. 

The restriction to the fundamental modular re- 
gion^U cancels the infrared singularity of the 
perturbation theory and makes the gluon to have 
the complex mass. The transverse gluon propa- 
gator in the continuum theory is defined as 



D MV {k) = - V e- lkx Tr{A^{x)A„{0)) 
n 

x— x,i 

= G(e)(s, v ~ k -^) . (i) 

where n is the dimension of the adjoint represen- 
tation of SU(N). Schematically 
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k 4 + K 4 2 y k 2 +iK 2 k 2 



(2) 



It was shown g that in the lattice simulation of 
the non-Abelian gauge theory, the Landau gauge 
fixing and the restriction of the gauge field yields 
the gluon propagator similar to that of Gribov. 
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2. The global gauge fixing and the gluon 
propagator 

We performed quenched (4 3 x 8 and 8 3 x 16) 
and unquenched (4 3 x 8) lattice QCD simulation 
in the Landau gauge whose algorithm is presented 
in HQ]. The number of configuration is 100 for the 
three cases. 

After the Landau gauge fixing with an accuracy 
Maxl^ylpl < 10~ 4 which is roughly equivalent 
to (d^A^) 2 < 10~ 10 , we fix the global gauge of 
the field A^{x) such that A = £ A a ^(x)X a of 
a sample is diagonalized by an SU(3) matrix g : 
g^Ag = diag(g 1 ,g 2 ,g 3 ), (gi > 92 > 9a). 

We obtain the new data sets A 9 ^(x) = 
g^A^(x)g and perform the Fourier transform 



,J2 x e- ik < x+e »WA 9a (x). 



Using the projection operator on a lattice in 
the momentum space 



Pnv(nk) = $, 



c p (n fc )c t ,(n fc ) 
c{n k ) 2 



(3) 



where c^(n k ) = 2sin^, c(n k ) 2 = Y^=i c n( n k) 2 
and fc M = ( nfc = 0, • • • , - 1) is the avail- 



able momentum on the lattice in the \x— direction, 
the gluon propagator is calculated from 



D(k) 



Y,P^k)((A 9 \(n k )A 9 \(n k y 



(A 9 %(n k ))(A 9 \(n k y)) 



(4) 



where (A ga fj,(n k )) is the average over samples of 
the globally gauge fixed field A^^rik). 
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We define the gluon propagator as a function 



of k = y'E^ K where fc M = Min[n fc , - n k }^ 

and not a function of k — ^/X^^sin-^-) 2 . 

In the case of 8 3 x 16 lattice there appear 
4096 independent four-momentum fc M . For a fixed 
k, the data of D(k) from < n k < and 
~2~ ^ "fe < which satisfies — n' k — are 
plotted at the same abscissa and in the medium 
momentum region the difference of the data at n' k 
and rik makes the plot of D(k) scattered. This 
effect could be attributed to the finite size effect 
and in larger lattice the effect becomes smaller IBJ. 

We adopt the data selection prescription sim- 
ilar to that of ||, i.e. select data such that 
\fiki — nkj\ < 1 where i and j run from 1 to 4, 
or choose the momentum k such that it is almost 
diagonal in the A 3 ^ lattice space. In fig. 1 we 
present the result of D (fc) . The data clearly shows 
the infrared suppression of the gluon propagator. 
The suppression remains even without the sub- 
traction of (A ga M (n k )) (A 9a v {n k )*) term. 

We performed the analysis also for k along the 
spacial axisQ. The results of 8 3 x 16 lattice sug- 
gests that the rotational symmetry is not so good. 
Although in the quenched 4 3 x 8 lattice the gluon 
propagator did not show the significant suppres- 
sion and remained finite, in the unquenched 4 3 x 8 
lattice similar suppression was observed |6j . In the 
SU(2) case similar effect is also observed^. 

3. Parametrization of the gluon propaga- 
tor 

In order to get insights on the complex mass of 
the confined gluon, we performed an analysis of 
the Fourier transform of the correlation function 

D T (k,t) = i-^ e 4k - x Tr(A 3 (x,r+t)^(0,T)) c (5) 

" jTx 

where the suffix c indicates that the connected 
part is taken after the global gauge fixing and 



— ik-s 



(6) 



In the Stingl's factorizing-denominator rational 
approximants (FDRA) method S, the transverse 
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Figure 1. The gluon propagator D(k) of quenched 
f3 = 5, 8 3 x 16 lattice. 



gluon propagator is expressed as 

D T (p2)[r=3] = P _j_ L 



p 2 + u + A 2 p 2 + t> + A 2 



+ c.c. (7) 



which is compared with its lattice Fourier trans- 
form 

GT(k, t) c = Re[(ao + zai)cosh((a 2 + ia 3 )(t — 4)) 

+ (04 + m 5 )cosh((a 6 + ia 7 )(t - 4))] 
pa 



Re 



L sinh(Mia)sinh(^Mia) 



-cosh(Mia(t - 4)) 



-cosh(M 2 a(* - 4))] (8) 



sinh(M 2 a)sinh(^M 2 a) 

where the gluon masses are defined as M\a = 
0-2 + *fl3, M2CL = a>6 + ia 7 and 

pa = (ao-Mai)sinh(a 2 +ia3)sinh(iV7'(a2-Ma3)/2),(9) 
ra = (a4+ia5)sinh(ae+iti7)sinh(iVT(ci6+ia7)/2),(10) 
u+A 2 = 4sinh 2 ((a 2 + io 3 )/2), (11) 



v + A 2 = 4sinh 2 ((a 6 + ia 7 )/2). 



(12) 



The gluon propagator D(p 2 ) derived in the 
FDRA method also shows infrared suppression. 
We observe that as the lattice size becomes large, 
the position of zero of D(p 2 ) approaches to p 2 = 
and that the peak near p 2 — becomes sharp. 
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Figure 2. The gluon propagator Gt(0,<) of 
quenched j3 = 5, 8 3 x 16 lattice. 
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Figure 3. The lattice Fourier transformed gluon 
propagator D(p 2 ). 



4. Discussion and outlook 

We performed the Langevin simulation of the 
quenched and unquenched lattice QCD in the 
Landau gauge, using the natural definition of the 
gauge field. 

The position of the zero of the Fourier trans- 
form of the connected part of the gluon propa- 
gator Dr(p 2 ) approaches p 2 = as the lattice 
size increases, which suggests the possibility of 
the realization of the confinement mechanism of 
the Gribov-Zwanziger theory. 

Our preliminary fitting results indicate that 
real part of the effective mass of the gluon is 
about 1.4a -1 for (3 = 4, 4 3 x 8 and about 0.75a _1 
for /3 = 5, 8 3 x 16. The renormalization and the 
anomalous dimension for the effective mass of the 
gluon is under investigation. 

This work was supported by High Energy Ac- 
celerator Research Organization, as KEK Super- 
computer Project(Project No. 98-34). 
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